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The short-lived 7^(892)* resonance provides an efficient tool to probe properties of the hot and 
dense medium produced in relativistic heavy-ion collisions. We report measurements of K* in ^/s^ 
= 200 GeV Au+Au and p+p collisions reconstructed via its hadronic decay channels /!'(892)**' — » /sTtt 
and 7^(892)** Klu^ using the STAR detector at RHIC. The K*° mass has been studied as a 
function of in mininmm bias p+p and central Au+Au collisions. The K* pr spectra for minimum 
bias p + p interactions and for Au+Au collisions in different centralities arc presented. The K* / K 
yield ratios for all centralities in Au+Au collisions are found to be significantly lower than the 
ratio in minimum bias p + p collisions, indicating the importance of hadronic interactions between 
chemical and kinetic freeze-outs. A significant non-zero K*° elliptic flow (^2) is observed in Au+Au 
collisions and compared to the K% and A V2- The nuclear modification factor of K* at intermediate 
Pt is similar to that of K%, but different from A. This establishes a baryon- meson effect over a mass 
effect in the particle production at intermediate pt (2 < pt < 4 GeV/c). 
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PACS numbers: 25.75.Dw, 25.75.-q, 13.75.Cs 



I. INTRODUCTION 

Lattice QCD calculations 1] predict a phase transition 
from hadronic matter to quark gluon plasma (QGP) at 
high temperatures and/or high densities. Matter under 
such extreme conditions can be studied in the labora- 
tory by colliding heavy nuclei at very high energies. The 
Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory provides collisions of heavy nuclei 
and protons at center of mass energies up to ^s„„ = 
200 GeV. The initial stage of these collisions can be de- 
scribed as the interpenetration of the nuclei with par- 
tonic interactions at high energy. With the interactions 
of the partons in the system, chemical and local ther- 
mal equilibrium of the system may be reached and the 
QGP may form. As the system expands and cools down, 
it will hadronize and chemically freeze-out. After a pe- 
riod of hadronic interactions, the system reaches the ki- 
netic freeze-out stage when all hadrons stop interact- 
ing After the kinetic freeze-out, particles free- 
stream towards the detectors where our measurements 
are performed. 

The typical lifetime of a resonance is a few fm/c, which 
is comparable to the expected lifetime of the hot and 
dense matter produced in heavy-ion collisions j^. In 
a hot and dense system, resonances are in close prox- 
imity with other strongly interacting hadrons. The in- 
medium effect related to the high density and/or high 
temperature of the medium can modify various reso- 
nance properties, such as masses, widths, and even the 
mass line shapes 0, S Hi- Thus, measurements of var- 
ious resonance properties can provide detailed informa- 
tion about the interaction dynamics in relativistic heavy- 
ion collisions ^ ^3 ■ Recent measurements [llj] by the 
FOCUS Collaboration for the K*'^ from charm decays 
show that the K*'^ mass line shape could be changed 
by the effects of interference from an s-wave and pos- 
sible other sources. Distortions of the line shape of 
have also been observed at RHIC in p + p and pe- 
ripheral Au-|-Au collisions p^a. Dynamical interactions 
with the surrounding matter 0, 1^ . Hp , interference be- 
tween various scattering channe ls llJ phase space dis- 
tortions 0, i, El IIIJS OMM and Bose- 
Einstein correlations |3, Mi llM l20| are possible 

explanations for the apparent modification of resonance 
properties. 

Resonance measurements in the presence of a dense 
medium can be significantly affected by two competing 
effects. Resonances that decay before kinetic freeze-out 
may not be reconstructed due to the rescattering of the 
daughter particles. In this case, the lost efficiency in the 
reconstruction of the parent resonance is relevant and de- 
pends on the time between chemical and kinetic freeze- 
outs, the source size, the resonance phase space distribu- 
tion, the resonance daughters' hadronic interaction cross- 



sections, etc. On the other hand, after chemical freeze- 
out, pseudo-elastic interactions |2ll | among hadrons in the 
medium may increase the resonance population. This 
resonance regeneration depends on the cross-section of 
the interacting hadrons in the medium. Thus, the study 
of resonances can provide an important probe of the time 
evolution of the source from chemical to kinetic freeze- 
outs and detailed information on hadronic interactions in 
the final stage. 

In this paper, we study the K{892)* vector meson with 
a lifetime of 4 fm/c. The kaon and pion daughters of the 
K* resonance in the hadronic decay channel K* — > Ktt 
can interact with other hadrons in the medium. Their 
rescattering effect is mainly determined by the pion-pion 
interaction total cross section p^ . which was measured 
to be significantly larger (factor ~5) than the kaon-pion 
interaction total cross section |2^. The kaon-pion in- 
teraction total cross section determines the regeneration 
effect that produces the K* resonance Thus, the 

final observable K* yields may decrease compared to the 
primordial yields, and a suppression of the K* / K yield 
ratio is expected in heavy- ion collisions. This K* yield 
decrease and the K* /K suppression compared to elemen- 
tary collisions, such as p-t-p, at similar collision energies 
can be used to roughly estimate the system time span 
between chemical and kinetic freeze-outs. Due to the 
rescattering of the daughter particles, the low px K* res- 
onances are less likely to escape the hadronic medium 
before decaying, compared to high px K* resonances. 
This could alter the K* transverse mass (tot) spectra 
compared to other particles with similar masses. 

The in-medium effects on the resonance production can 
be manifested in other observables as well. In a quark co- 
alescence scenario, the elliptic flow (W2), for non-central 
Au-f Au collisions, of the K* resonances produced at 
chemical freeze-out might be similar to that of kaons 2^. 
However, at low pt, the K* V2 may be modified by the 
rescattering effect discussed previously. This rescatter- 
ing effect also depends on the hadron distributions in the 
coordinate space in the system at the final stage. Thus, 
a measurement of the K* V2 at pt < 2 GeV/c compared 
to the kaon V2 may provide information on the shape of 
the fireball in the coordinate space at late stages. 

A study of the relation of the particle production to its 
intrinsic properties may reveal its production mechanism. 
The nuclear modification factor and V2 have been ob- 
served to be different between tt, K and p, A [Hl^^l. In 
a hydrodynamic limit, the transverse momentum spectra 
of produced particles are only determined by the velocity 
field and therefore the mass of the produced particle. In 
a quark coalescence model, particle production is related 
to its quark content. Since stable mesons (tt, K) are usu- 
ally lighter than stable baryons {p, A), the particle type 
is coupled with the mass. Detailed studies of K* (and/or 
(j)) can be of special importance, as its mass is close to 
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the mass of baryons {p, A) but it is a vector meson. In 
the intermediate px range 2 < pt < G GeV/c, identi- 
fied hadron V2 measurements have shown that the hadron 
V2 follows a simple scaling of the number of constituent 
quarks in the hadrons: V2{pt) = ™2{PT/'n), where n is 
the number of constituent quarks of the hadron and 
is the common elliptic flow for sing le quarks f27l|. There- 
fore, the V2 for the K* produced at hadronization should 
follow the scaling law with n — 2. However, for the K* 
regenerated through K* in the hadronic stage, 

V2 should follow the scaling law with n = 4 The 
measured K* V2 in the intermediate px region may pro- 
vide information on the K* production mechanism in 
the hadronic phase and reveal the particle production 
dynamics in general. It is inconclusive whether the dif- 
ference in the nuclear modification factor between K and 
A is due to a baryon-meson effect or simply a mass ef- 
fect ■ We can use the unique properties of the K* 
to distinguish whether the nuclear modification factor 
Raa or i?r7 p, defined later in the text (Section V.F.) and 
in [29j | and |27l | respectively, in the intermediate pt region 
depends on mass or particle species (i.e. meson/baryon). 
Specifically, we can compare the Rcp of K, K*, and 
A, which contain one strange valence quark and are in 
groups of {K, K*) and A as mesons vs baryon, or in 
groups of K and (if*. A) as different masses. 



duction, the events from minimum bias Au-|-Au collisions 
were divided into four centrality bins from the most cen- 
tral to the most peripheral coUisions: 0-10%, 10-30%, 30- 
50% and 50-80%, according to the fraction of the charged 
hadron reference multiplicity (defined in |31j) distribu- 
tion in all events. 




^ p (GeV/c) 



FIG. 1: (Color online) dE/dx for negative particles vs. mo- 
mentum measured by the TPC in Au-|-Au collisions. The 
curves are the Bethe-Bloch parametrization [3^ for different 
particle species. 



II. EXPERIMENT 

The data used in this analysis were taken in the second 
RHIC run (2001-2002) using the Solenoidal Tracker at 
RHIC (STAR) with Au-|-Au and p-l-p collisions at 
= 200 GeV. The primary tracking device of the STAR 
detector is the time projection chamber (TPC) which is a 
4.2 meter long cylinder covering a pseudo-rapidity range 
I77I < 1.8 for tracking with complete azimuthal coverage 
(A0 = 27r) 

In Au-|-Au collisions, a minimum bias trigger was de- 
fined by requiring coincidences between two zero degree 
calorimeters which are located in the beam directions at 
<2 mrad and measure the spectator neutrons. A cen- 
tral trigger corresponding to the top 10% of the inelastic 
hadronic Au-|-Au cross-section was defined using both the 
zero degree calorimeters and the scintillating central trig- 
ger barrel, which surrounds the outer cylinder of the TPC 
and triggers on charged particles in the midpseudorapid- 
ity (I77I < 0.5) region. In p -f p collisions, the minimum 
bias trigger was defined using coincidences between two 
beam-beam counters that measure the charged particle 
multiplicity in forward pseudorapidities (3.3 < \ri\ < 5.0). 

Only events with the primary vertex within ±50 cm 
from the center of the TPC along the beam line were se- 
lected to insure uniform acceptance in the 77 range stud- 
ied. As a result, about 2 x 10^ top 10% central Au-|-Au, 
2 X 10® minimum bias Au+Au, and 6 x 10® minimum 
bias p + p collision events were used in this analysis. In 
order to study the centrality dependence of the K* pro- 



In addition to momentum information, the TPC pro- 
vides particle identification for charged particles by mea- 
suring their ionization energy loss {dE/dx). The TPC 
measurement of dE/dx as a function of the momentum 
[p) is shown in Fig. ^ Different bands seen in Fig. 
represent Bethe-Bloch distributions [s^ folded with the 
experimental resolutions and correspond to different par- 
ticle species. Charged pious and kaons can be identified 
with their momenta up to about 0.75 GeV/c while pro- 
tons and anti-protons can be identified with momenta up 
to about 1.1 GeV/c. In order to quantitatively describe 
the particle identification, the variable N^^^ (e.g. pions) 
was defined as 



1 }q /^-^measured 

'R {dE/dx)^ '■ 



(1) 



where dE / dxmcasuvcd is the measured energy loss for a 
track, {dE/dx)^^ is the expected mean energy loss for 
charged pions with a given momentum, and R is the 
dE/dx resolution which varies between 6% and 10% from 
p + p to central Au+Au events and depends on the char- 
acteristics of each track, such as the number of dE/dx 
hits for a track measured in the TPC, the pseudorapid- 
ity of a track, etc. We construct N^k in a similar way 
for the charged kaon identification. Specific analysis cuts 
(described later) were then applied on N^-tt E^nd N^-k in 
order to quantitatively select the charged pion and kaon 
candidate tracks. 
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III. PARTICLE SELECTIONS 



In this analysis, 
K{892)*° ^ K+ 



the hadronic decay channels of 
TT-, K{892)*o K-TT+ and 

* were measured. In the following, 
and the term K* 



or K 



*o 



is:(892)*± ^ K°Tr 
the term K*° s tands for 
stands for K*^, K*'^ or K*^, unless otherwise specified. 

Since the K* decays in such short time that the daugh- 
ters seem to originate from the interaction point, only 
charged kaon and charged pion candidates whose dis- 
tance of closest approach to the primary interaction ver- 
tex was less than 3 cm were selected. Such candidate 
tracks are defined as "primary tracks" . The charged K* 
first undergoes a strong decay to produce a Kg and a 
charged pion herein labeled as the K*^ daughter pion. 
Then, the produced Kg decays weakly into tt+tt" with 
CT = 2.67 cm. Two oppositely charged pions from the Kg 
decay are called as the K*"^ granddaughter pions. The 
charged daughter pion candidates were selected from pri- 
mary track samples and the A'^ candidates were selected 
through their decay topology. 

In Au-|-Au collisions, charged kaon candidates were 
selected by requiring \NaK\ < 2 while a looser cut 
|7Vo-7r| < 3 was applied to select the charged pion can- 
didates to maximize the statistics for the K*'^ analysis. 
Such Ncr cuts can only ambiguously select the kaons and 
pions if applied to the tracks with their momenta beyond 
the momentum range specified earlier. However, these 
cuts help to significantly reduce the background. In or- 
der to avoid the acceptance drop in the high rj range, all 
kaon and pion candidates were required to have I77I < 0.8. 
Kaon and pion candidates were also required to have at 
least 15 fit points (number of measured TPC hits used 
in track fit, maximum 45 fit points) to assure track fit- 
ting quality and good dE/dx resolution. For all the track 
candidates, the ratio between the number of TPC track 
fit points over the maximum possible points was required 
to be greater than 0.55 to avoid selecting split tracks. To 
maintain reasonable momentum resolution, only tracks 
with pt larger than 0.2 GeV/c were selected. 

In p + p collisions, enough data were available to pre- 
cisely measure the K*^ mass, width, and invariant yield 
as a function of pt- As statistics was not an issue for this 
analysis, only kaon candidates with p < 0.7 GeV/c were 
used to ensure clean identification. This kaon momen- 
tum cut helped minimize contamination from misidenti- 
fied correlated pairs and thus reduce the systematic un- 
certainty. In the case of the pion candidates, the same 
p and Pt cuts as used in Au-|-Au collisions were applied. 
Charged kaon and pion candidates were selected by re- 
quiring |A^cr7r,i<-| < 2 to reduce the residual background. 
All other track cuts for both kaon and pion candidates 
were the same as for Au-|-Au data. 

The K*^ was measured only in minimum bias p + p 
interactions and in peripheral 50-80% Au-|-Au collisions. 
Daughter pions for the K*^ reconstruction were required 
to originate from the interaction point and pass the 
same cuts as used for the K*'^ analysis in p + p colli- 




ji+r Inv. Mass (GeV/c ) 



FIG. 2; Kg signal observed in the tt+tt" invariant mass dis- 
tribution reconstructed from the decay topology method via 
Kg — > n^Tv~ in p + p collisions. The dashed curve depicts the 
Gaussian fit function plus a linear function representing the 
background. 



sions. The A"" was reconstructed by the decay topology 
method [33, yj] . The granddaughter charged pion candi- 
dates were selected from global tracks (tracks do not nec- 
essarily originate from the primary collision vertex) with 
a distance of closest approach to the interaction point 
greater than 0.5 cm. Candidates for the granddaugh- 
ter charged pions were also required to have at least 15 
hit points in the TPC with p > 0.2 GeV/c. Oppositely 
charged candidates were then paired to form neutral de- 
cay vertices. The distance of closest approach for each 
pair was required to be less than 1.0 cm and the neu- 
tral decay vertices were required to be at least 2.0 cm 
away from the primary vertex to reduce the combinato- 
rial background. The reconstructed Kg momentum vec- 
tor was required to point back to the primary interac- 
tion point within 1.0 cm. Only the Kg candidates with 
TT+TT^ invariant mass between 0.48 and 0.51 GeV/c^ were 
selected. When the Kg candidate was paired with the 
daughter pion to reconstruct the charged K* , tracks were 
checked to avoid double-counting among the three tracks 
used. Fig. 13 shows the Kg signal observed in the tt+tt" 
invariant mass distribution in p + p collisions. The Gaus- 
sian width of the above Kg signal is around 7 MeV/c^ 
which is mainly determined by the momentum resolu- 
tion of the detector. Due to detector effects, such as the 
daughter tracks' energy loss in the TPC, etc., the Kg 
mass is shifted by —3 MeV/c^. The measured Kg mass 
and width agree well with Monte Carlo (MC) simula- 
tions, which included the finite momentum resolution of 
the detector and the daughter tracks' energy loss in the 
TPC. 



The Ktt pairs with their parent rapidity (y) of \y\ < 0.5 
were selected. All the cuts used in this K* analysis are 
summarized in Table After all the above mentioned 
cuts have been applied, the K* reconstruction efficiencies 
multiplied by the detector acceptance are shown in Fig.|51 
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TABLE I: List of track cuts for charged kaon and charged pion and topological cuts for neutral kaon used in the K* analysis 
in Au+Au and p + p collisions. decayLength is the decay length, dcaDaughters is the distance of closest approach between 
the daughters, dcaVOPrmVx is the distance of closest approach between the reconstructed Kg momentum vector and the 
primary interaction vertex, dcaPosPrmVx is the distance of closest approach between the positively charged granddaughter 
and the primary vertex, dcaNegPrmVx is the distance of closest approach between the negatively charged granddaughter and 
the primary vertex, M^^o is the K% invariant mass in GeV /(? , N FitPnts is the number of fit points of a track in the TPC, 
NTpcHits is the number of hits of a track in the TPC, MaxPnts is the number of maximum possible points of a track in the 
TPC, and DC A is the distance of closest approach to the primary interaction point. 



Cuts 


Au+Au 


p + p 


Daughter tt* 




N„K 


(-2.0, 2.0) 


(-2.0, 2.0) 




decayLength > 2.0 cm 




(-3.0, 3.0) 


(-2.0, 2.0) 


(-2.0, 2.0) 


dcaDaughters < 1.0 cm 


Kaon p (GeV/c) 


(0.2, 10.0) 


(0.2, 0.7) 




dcaVOPrmVx < 1.0 cm 


Kaon pt (GeV/c) 


(0.2, 10.0) 


(0.2, 0.7) 




dcaPosPrmVx > 0.5 cm 


Pion p (GeV/c) 


(0.2, 10.0) 


(0.2, 10.0) 


(0.2, 10.0) 


dcaN egPrmV x > 0.5 cm 


Pion PT (GeV/c) 


(0.2, 10.0) 


(0.2, 10.0) 


(0.2, 10.0) 


M^o (GeV/c^): (0.48, 0.51) 


NFitPnts 


> 15 


> 15 


> 15 


7r+: NTpcHits > 15 


N FitPnts /M ax Puts 


> 0.55 


> 0.55 


> 0.55 


TT": NTpcHits > 15 


Kaon and pion rj 


\r]\ < 0.8 


\V\ < 0.8 


|r/l < 0.8 


TV+: p> 0.2 GeV/c 


DC A (cm) 


< 3.0 


< 3.0 


< 3.0 


p> 0.2 GeV/c 


Pair (Kn) y 






\y\ < 0.5 





IV. EXTRACTION OF THE A* SIGNAL 

In Au+Au collisions, up to several thousand charged 
tracks per event originate from the primary collision 
vertex. The daughters from K* decays are topologi- 
cally indistinguishable from other primary particles. The 
measurement was performed by calculating the invariant 
mass for each Ktt pair in an event. The K^n^ invari- 
ant mass distribution is shown in Fig. |3| as open circles. 
The unlike-sign Kit invariant mass distribution derived 
in this manner was mostly from random Ktt combina- 
torial pairs. The signal to background is between 1/200 
for minimum bias Au+Au and 1/10 for minimum bias 
p + p. The overwhelming combinatorial background dis- 
tribution can be obtained and subtracted from the unlike- 
sign Kt: invariant mass distribution in two ways: 

• the mixed-event technique: reference background 
distribution is built with uncorrelated unlike-sign 
kaons and pions from different events; 

• the like-sign technique: reference background dis- 
tribution is made from like-sign kaons and pions in 
the same event. 

The mixed-event technique has been successfully used 
in the measurement of resonances at RHIC, such as the 
K{892)*" in Au+Au coUisions at = 130 GeV 

and the (j) in Au+Au collisions at y/s^ = 130 and 
200 GeV 113' This technique was also used in 
tlie measurement of A production in Au+Au collisions 
^/s^ = 130 GeV, and the result s agree well with 
those from the decay topology method |34Ll38| . The like- 
sign technique has been successfully applied in measur- 
ing p(770)° — > TT+TT" production in p + p and peripheral 
Au+Au coUisions at = 200 GeV at RHIC (IJJ. 



A. Mixed-Event Technique 

In order to subtract the uncorrelated pairs from the 
unlike-sign Kn invariant mass distribution obtained from 
the same events, an unlike-sign Kir invariant mass spec- 
trum from mixed events was obtained. In order to keep 
the event characteristics as similar as possible among dif- 
ferent events, the whole data sample was divided into 10 
bins in charged particle multiplicity and 10 bins in the 
collision vertex position along the beam direction. Only 
pairs from events in the same multiplicity and vertex po- 
sition bins were selected. 




0.6 0.8 1 1.2 1.4 



Ktt Inv. Mass (GeV/c^) 

FIG. 3: The unlike-sign Kn invariant mass distribution (open 
symbols) and the mixed-event Kn invariant mass distribution 
after normalization (solid curve) from minimum bias Au+Au 
collisions. 

In the unlike-sign invariant mass distribution from an 
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event, K^tt^ and K^n^ pairs were sampled which in- 
clude the desired K* signal and the background. In 
the mixed-event spectrum, K^n^ , K^n^, K^tt^ , and 
K~7r^ pairs were sampled for the background estimation. 
The subscripts 1 and i correspond to event numbers with 
I 7^1. The number of events to be mixed was chosen to 
be 5, so that the total number of entries in the mixed- 
event invariant mass distribution was ~10 times that of 
the total number of entries in the distribution from the 
same events. Thus the mixed-event spectrum needs to 
be normalized in order to subtract the background in the 
unlike-sign spectrum. Since the Kit pairs with invariant 
mass greater than 1.1 GeV/c^ are less likely to be cor- 
related in the unlike-sign distribution, the normalization 
factor was calculated by taking the ratio between the 
number of entries in the unlike-sign and the mixed-event 
distributions for invariant mass greater than 1.1 GeV/c^. 
The solid curve in Fig. |3| corresponds to the mixed-event 
Kn pair invariant mass distribution after normalization. 
The mixed-event distribution was then subtracted from 
the unlike-sign distribution as follows: 

NK^oim) = N^+^- (m) -|- iV^-^+ (m) 

6 

+Nj,+^-{m)+Nj,-^+{m)], (2) 

where N is the number of entries in a bin with its center 
at the Kt: pair invariant mass m and R is the normaliza- 
tion factor. After the mixed-event background subtrac- 
tion, the K*^ signal is visible as depicted by the open 
star symbols in Fig. |S1 

B. Like-Sign Technique 

The like-sign technique is another approach to sub- 
tract the background of non-correlated pairs from the 
unlike-sign K-k invariant mass distribution from the same 
events. The uncorrelated background in the unlike-sign 
Kit distribution was described by using the invariant 
mass distributions obtained from uncorrelated X+7r+ 
and K~Tr~ pairs from the same events. 

In the unlike-sign Ktt invariant mass spectrum, K^tt^ 
and K^TT^ pairs were sampled. K^n^ and K^n^ pairs 
were sampled in the like-sign Ktt invariant mass distribu- 
tion. Since the number of positive and negative particles 
may not be the same in relativistic heavy-ion collisions, 
in order to correctly subtract the subset of non-correlated 
pairs in the unlike-sign Kt: distribution, the like-sign Kn 
invariant mass distribution was calculated as follows: 

iVLiko-Sign(™) = 2 X ^Nj^+^+{m) X N^-^-{m), (3) 

where N is the number of entries in a bin with its center 
at the Kt: pair invariant mass m. The unlike-sign and 
the like-sign invariant mass distributions are shown in 



Fig. 0] The like-sign spectrum was then subtracted from 
the unlike-sign distribution: 

NK'o{m) = Nj^+^- (to) -I- Nj^-^+ (m) - Nukc-Sign{m) . 

. . (4) 

The like-sign background subtracted Kn invariant mass 
distribution corresponds to the solid square symbols in 
Fig. [Sj where the K*'^ signal is now visible. 

Compared to the mixed-event technique, the like-sign 
technique has the advantage that the unlike-sign and like- 
sign pairs are taken from the same event, so there is no 
event structure difference between the two distributions 
due to effects such as elliptic flow. The short-coming 
of this technique is that the like-sign distribution has 
larger statistical uncertainties compared to the mixed- 
event spectrum, since the statistics in the mixed-event 
and like-sign techniques are driven by the number of 
events mixed and the number of kaons and pions pro- 
duced per event, respectively Therefore, in this 
analysis, the mixed-event technique was used to recon- 
struct the K* signal whereas the like-sign technique was 
used to study the sources of the residual background un- 
der the K*'^ peak after mixed-event background subtrac- 
tion as discussed in details in the following text. 




0.6 0.8 1 1.2 1.4 



Ktt Inv. Mass (GeV/c^) 

FIG. 4: The unlike-sign Ktt invariant mass distribution (open 
symbols) and the like-sign Ktt invariant mass distribution 
(solid curve) from minimum bias Au-|-Au collisions. 



C. Describing the Residual Background 

The unlike-sign Ktt invariant mass distribution after 
mixed-event background subtraction is represented by 
the open star symbols in Fig. where the K*^ signal 
is clearly observed. The mixed-event technique removes 
only the uncorrelated background pairs in the unlike-sign 
spectrum. As a consequence, residual correlations near 
the K*^ mass range were not subtracted by the mixed- 
event spectrum. This residual background may come 
from three dominant sources: 
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• elliptic flow in non-central Au+Au collisions; 

• correlated real Kt: pairs; 

• correlated but misidentified pairs. 

The overlapping region of non-central Au-f Au colli- 
sions has an elliptic shape in the plane perpendicular 
to the beam axis. Each non-central Au-f Au event has 
a unique reaction plane angle. The azimuthal distribu- 
tions for kaons and pions may be different for different 
events. Thus, the unlike-sign Kn pair invariant mass 
spectrum may have a different structure than the mixed- 
event invariant mass distribution. This structural dif- 
ference may lead to a significant residual background in 
the unlike-sign Kn invariant mass spectrum after mixed- 
event background subtraction [40l |. 

In the like-sign technique, the unlike-sign Kn spec- 
trum and the like-sign distribution are obtained from 
the same events. Therefore, no correlations due to ellip- 
tic flow should be present in the unlike-sign Kn invari- 
ant mass spectrum after like-sign background subtrac- 
tion. In Fig. O the solid square symbols represent the 
unlike-sign Kt: invariant mass distribution after like-sign 
background subtraction. The amplitude of the residual 
background below the peak after the like-sign background 
subtraction is about a factor of 2 smaller than after the 
mixed-event background subtraction, while the ampli- 
tude of the K*'^ signal remains the same. This indicates 
that part of the residual background in the spectrum af- 
ter mixed-event background subtraction was induced by 
elliptic azimuthal anisotropy. 
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FIG. 5: The Kn invariant mass distributions after event- 
mixing background subtraction (open star symbols) and like- 
sign background subtraction with different daughter momen- 
tum cuts (0.2 < Kaon and Pion p < 10 GeV/c for filled square 
symbols, 0.2 < Kaon p < 0.7 GeV/c and 0.2 < Pion p < 10 
GeV/c for open triangle symbols) demonstrating the sources 
of the residual background in minimum bias Au-|-Au colli- 
sions. The open triangle symbols have been scaled up by a 
factor of 3 in order to increase the visibility. The arrow de- 
picts the standard K'° mass of 896.1 MeV/c^ 



In the K*'^ analysis in Au-t-Au collisions, since the 
kaons and pions are selected with 0.2 < p < 10.0 GeV/c, 
a pion (kaon) with p > 0.75 GeV/c may be misidenti- 
fied as a kaon (pion). A proton with p > 1.1 GeV/c 
may be misidentified as either a kaon or a pion, or both, 
depending on whether kaons or pions are being selected. 
Electrons and positrons which cross the kaon (pion) band 
in the dE/dx plot shown in Fig. ^^lay be misidentifed 
as kaons (pions). Thus, the daughters from — > tt+tt", 
(j) K~^K~, A — > n~p, etc. could be falsely identified as 
a Kn pair if the daughter momenta are beyond the par- 
ticle identification range. The invariant mass calculated 
from these misidentified pairs cannot be subtracted away 
by the mixed-event background and remains as part of 
the residual background. 

In Fig. El the open triangle symbols correspond to the 
unlike-sign Kn spectrum after like-sign background sub- 
traction with 0.2 < p < 0.7 GeV/c and 0.2 < p < 10.0 
GeV/c for the kaon and the pion, respectively. These 
momentum cuts allow only correlated Kn real pairs and 
pairs in which a kaon or a proton was misidentified as 
a pion to contribute to the background subtracted spec- 
trum. Compared to the solid square symbols in Fig. \S\ 
the residual background represented by the open trian- 
gle symbols is reduced by a factor of 6 and the K*° sig- 
nal is a factor of 2 smaller. This indicates that particle 
misidentification of the decay products of p, co, r], Kg, 
A, etc. indeed causes false correlations to appear in the 
background subtracted distribution. 

Correlated real Kn pairs from real particle decays, 
such as higher mass resonant states in the K — n sys- 
tem and particle decay modes with three or more daugh- 
ters where two of them are a Kn pair, as well as the 
nonresonant K — n s-wave correlation also contribute to 
the unlike-sign Kn spectrum. These correlated Kn pairs 
contribute to the residual background, since they are not 
present in the like-sign and mixed-event distributions. 
There is no efficient cut to remove these real correlations 
from the residual background. 



V. ANALYSIS AND RESULTS 

A. K*° Mass and Width 

Figure El depicts the mixed-event background sub- 
tracted Kn invariant mass distributions {Mk-k) inte- 
grated over the K* pT for central Au-|-Au (upper panel) 
and for minimum bias p + p (lower panel) interactions. 
The mass of the K*'^ was fit to the function: 



BW X PS + RBG, 



(5) 



where BW is the relativistic p-wave Breit-Wigner func- 
tion [H: 



BW 



Mk^TMo 



(6) 
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PS is the Boltzmann factor 

Mktt 



PS = 



X exp 



mm 

Tfo 



(7) 



that accounts for phase space, and RBG is the hnear 
function: 



RBG = 



bM. 



Kit- 



(8) 



that represents the residual background. Within this 
parametrization, Ttn is the temperature at which the res- 
onance is emitted Ig and 



[Ml^-Ml-Mlf 



{Mi 



M2 - M|.)2 



TT K 



3/2 



(9) 



is the momentum dependent width j41| . In addition, Mq 
is the K* mass, Fq is the K* width, pr is the K* trans- 
verse momentum, M^^ is the pion mass, and Mj^ is the 
kaon mass. 
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FIG. 6: The Ktt invariant mass distribution integrated over 
the K* pt for central Au+Au (upper panel) and minimum 
bias p + p (lower panel) interactions after the mixed-event 
background subtraction. The solid curves are the fits to Eq.|3 
with Tfo = 120 MeV and pr = 1.8 GeV/c for central Au+Au 
and Tfo = 160 MeV and pr = 0.8 GeV/c for p + p, respec- 
tively. The dashed lines are the linear function representing 
the residual background. 

The PS factor accounts for K* produced through kaon 
and pion scattering, or K+tt K* K+n. In Au-|-Au 
collisions, the thermal freeze-out temperature Tfo = 90 
MeV was measured at STAR l42|. However, resonances 



can be produced over a range of temperature inside the 
hadronic system and not all resonances are emitted at the 
point where the system freezes out at Tfo = 90 MeV. As 
a result, the temperature chosen in the PS factor was 120 
MeV according to The temperature of Tfo = 90 MeV 
was also used to estimate the systematic uncertainties 
which are about 1.5 MeV/c^ for masses and 5 MeV/c^ 
due the choice of Tfo- In p+p collisions, particle produc- 
tion is well reproduced by the statistical model |43] with 
Tfo = 160 MeV and therefore this was the temperature 
used in the PS factor, pr = 1.8 GeV/c and 0.8 GeV/c 
were chosen in the PS factor for the Au-|-Au and p + p 
collisions respectively which are the centers of the entire 
measured pr ranges (0.4 < px < 3.2 GeV/c for Au+Au 
and Pt <1.6 GeV/c for p+p). 
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FIG. 7: The K*'^ mass (upper panel) and width (lower panel) 
as a function of pr for minimum bias p + p interactions 
and for central Au+Au collisions. The solid straight lines 
are the standard K'^ mass (896.1 MeV/c^) and width (50.7 
MeV/c^) |3^. respectively. The dashed and dotted curves are 
the MC results in minimum bias p+p and for central Au+Au 
collisions, respectively, after considering detector effects and 
kinematic cuts. The grey shadows (caps) indicate the system- 
atic uncertainties for the measurement in minimum bias p + p 
interactions (central Au+Au collisions). 

Mixed-event background subtracted Kt: invariant 
mass distributions were obtained for different pr bins, 
and each pr bin was fit to Eq. [S] with the if*" mass, 
width, and uncorrected yield as free parameters. The 
X^/ndf of the fit varies between 0.6 and 1.7 for all pT 
bins except for two pr bins (3.8 for the 2.0 < pr < 2.4 
GeV/c bin and 2.6 for the 2.4 < pr < 2.8 GeV/c bin) 
in the central Au+Au data, where the uncertainties of 
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the mass and width values are not well constrained. Fig- 
ureElshows the K*'^ mass (upper panel) and width (lower 
panel) for central Au+Au and for minimum bias p + p 
interactions as a function of the K*'^ pT- MC calcula- 
tions for the K*^ mass and width were obtained by sim- 
ulating if *° with standard mass and width values [s^ ] 
and passing them through the same reconstruction steps 
and kinematic cuts as the real data. The results from 
such simulations are also depicted in Fig. [3 The devi- 
ations between the MC results and the standard mass 
and width values are mainly due to the kinematic cuts 
(track p and pr cuts, etc.). For example, the px > 0.2 
GeV/c cut results in the rise of the mass at low px and 
the kaon p < 0.7 GeV/c cut mp+p causes the rise of the 
mass and the drop of the width at higher px- Our MC 
studies indicate that the deviations induced by kinematic 
cuts are not sufficient to explain the mass shift seen in 
the data. 

The systematic uncertainties in the K*^ mass and 
width for the measurement in minimum bias p + p in- 
teractions were evaluated bin-by-bin by varying the par- 



ticle types (either K 



*o 



K*^ ), the methods in the 



background subtraction (mixed-event or like-sign), the 
residual background functions (exponential or second or- 
der polynomial functions), the dynamical cuts, the track 
types (primary tracks or global tracks), and by consid- 
ering the detector effects (different TPC magnetic field 
directions, different sides of the TPC detector, etc.). Due 
to the limited statistics, the systematic uncertainties (3.1 
MeV/c^ for masses and 14.9 MeV/c^ for widths) in cen- 
tral Au-fAu interactions were only estimated using the 
entire measured range (0.4 < px < 3.2 GeV/c) follow- 
ing the above steps. More detailed discussions about the 
systematic uncertainty studies can be found in [s^- In 
minimum bias p-|-p interactions, the K*^ masses at low 
Px (first 2-3 data points) are lower than the MC results 
at a 2-3(7 level. 



B. rriT and px Spectra 

Mixed-event background subtracted Ktt invariant 
mass distributions were obtained for different px bins, 
and each px bin was fit to function: 



SBW + RBG, 



(10) 



where SBW is the non-relativistic Breit-Wigner function 



SBW = 



Fn 



(11) 



(Mk^ - Mo)2 + F2/4 

and RBG is the linear function from Eq. |H1 that rep- 
resents the residual background. The fit sensitivity to 
statistical fluctuations in the K* raw yield was reduced 
by fixing the mass and width in the fit according to the 
values obtained from the free parameter fit with the same 
simplified BW function. The x^/^df of the fit varies be- 
tween 0.7 and 1.8 for all px bins except for two px bins 



(~3.0 for the 2.0 < px < 2.4 GeV/c bin and ~2.7 for the 
2.4 <px < 2.8 GeV/c bin) in Au-hAu data. The K* raw 
yield was also obtained by fitting the data to the BW 
function from Eq. with all parameters free in the fit. 
The difference in the raw yields between the two fit func- 
tions was included in the systematic uncertainties. The 
KgTT^ invariant mass distribution fit to Eq. ^| after the 
mixed-event background subtraction is shown in Fig. |S1 
for minimum bias p + p collisions (upper panel) and for 
the 50-80% of the inelastic hadronic Au-|-Au cross-section 
(lower panel). 




K"7t± Inv. Mass (GeV/c ) 



FIG. 8: The K^giv^ invariant mass distribution integrated over 
the K*^ pt for minimum bias p + p collisions (upper panel) 
and for the 50-80% of the inelastic hadronic Au+Au cross- 
section (lower panel) after the mixed-event background sub- 
traction. The solid curves are fits to Eq. 1101 and the dashed 
lines are the linear function representing the residual back- 
ground. 

About 6xl0^ 2x10^ and 5.6x10^ K*" signals were 
reconstructed from top 10% central Au+Au, minimum 
bias Au+Au and minimum bias p + p collisions respec- 
tively while about 1.2x10^ and 10^ K*^ were observed 
in the 50-80% Au+Au and minimum bias p + p collisions 
respectively. The A'*" and K*^ raw yields obtained for 
different px bins in Au+Au and minimum bias p+p colli- 
sions were then corrected for the detector acceptance and 
efficiency (shown in Fig. I^J determined from a detailed 
simulation of the TPC response using GEANT 01 • The 
corresponding branching ratios were also taken into ac- 
count. In addition, the yields in p + p were corrected for 
the collision vertex finding efficiency of 86%. 

The transverse mass (my) distributions of the midra- 
pidity {K*'^ + K*'^)/2 invariant yields in central Au+Au, 
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FIG. 9: The -K"*" and K*^ reconstruction efficiency multiplied 
by the detector acceptance as a function of pr in minimum 
bias p + p and different centralities in Au+Au collisions. 
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four different centralities in minimum bias Au+Au, and 
minimum bias p+p collisions are depicted in Fig. llUI The 
(A'*^ + K*^)/2 invariant yields for the most peripheral 
50-80% Au+Au collisions are also shown for comparison. 
The K*° invariant yield [(P N / {2TTmTdydniT)] distribu- 
tions were fit to an exponential function: 



1 



dN 



1 



2iTmT dydrriT 



exp 



dy 27rr(mo + T) 
[rnr - mo) 



T 



(12) 



where dN/dy is the K*'^ yield at \y\ < 0.5 and T is the 
inverse slope parameter. The extracted dN/dy and T 
parameters are listed in Tabled The systematic uncer- 
tainties on the K*'^ dN/ dy and T in Au+Au and p + p 
collisions were estimated by comparing different Breit- 
Wigner functions, particle types (either K*'^ or K*^), 
residual background functions (exponential or second or- 
der polynomial functions), dynamical cuts, and by con- 
sidering the detector effects. More detailed discussions 
about the systematic uncertainty studies can be found 
in [39| . The K**^ invariant yield increases from p + p 
collisions to peripheral Au+Au and to central Au+Au 
collisions. The inverse slope of the if*" spectra for all 
centrality bins of Au+Au collisions is significantly larger 
than in minimum bias p + p collisions. 

Theoretical calculations [43 indicate that in p + p col- 
lisions, particle production is dominated by hard pro- 
cesses for pt above 1.5 GeV/c while soft processes dom- 
inate at low Pt- Thus in the K* px spectrum, a power- 
law shape for pT above 1.5 GeV/c and an exponential 
shape at lower pT should be expected. In minimum bias 
p+p collisions, due to the cut on the kaon daughter 
of p < 0.7 GeV/c, only the K*^ spectrum for pr < 1.6 
GeV/c was measured. As a result, the K*'^ rriT spectrum 
in minimum bias p + p collisions can be well described by 
the commonly used exponential function, as shown in 
Fig. ^1 The K* pr spectrum can be extended to higher 



FIG. 10: The {K* + K*)/2 invariant yields as a function of 
rriT — mo for \y\ < 0.5 from minimum bias p + p and different 
centralities in Au+Au collisions. The top 10% central data 
have been multiplied by two for clarity. The lines are fits to 
Eq. 1121 The errors shown are the quadratic sum of statistical 
and systematic (in the level of 10%) uncertainties. 



TABLE II: The K*° dN/dy and T for |y| < 0.5 from central 
Au+Au, four different centralities in minimum bias Au+Au, 
and minimum bias p+p collisions. The first error is statistical, 
the second is systematic. 



top 10% central 
0-10% 
10-30% 
30-50% 
50-80% 
p + p 



dN/dy 



T (MeV) 



10.18+0.46+1.88 
10.48+1.45+1.94 
5.86+0.56+1.08 
2.81+0.25+0.52 
0.82+0.06+0.15 
(5.08+0.17+0.61) X 10" 



427+10+46 
428+31+47 
446+23+49 
427+18+46 
402+14+44 
223+8+9 



Pt by measu ring the K*^ signals. Figure ITTI shows the 
(A'*" +^)/2 and {K*+ + K*-)/2 invariant yields for 
\y\ < 0.5 as a function of pr- The dotted curve in this 
figure is the fit to the power-law function: 



1 



d'^N _ dN2{n- l)("-2) 



2TrpT dydpT 



1 



dy Tr{n - Syipr) 
Pt 

(pT)(n-3)/2 



(13) 



where n is the order of the power law and (px) is the 
average transverse momentum. The data were fit for 
Pt > 0.5 GeV/c. The power-law fit does not reproduce 
the two first pr bins (0.0 < pr < 0.2 GeV/c and 0.2 
< Pt < 0.4 GeV/c) since at low px particle production 
may be dominated by soft processes. From the power-law 
fit, the /ndf is 0.93. The dashed curve in Fig.^Jis the 
K*^ spectrum fit to the exponential function from Ea. ll2l 
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and then extrapolated to higher px- The data could not 
be described by this exponential fit indicating that hard 
processes dominate the particle production for pt > 1.5 
GeV/c. Some model |4g suggests to use the Levy func- 
tion in Equation 1141 to represent the pr spectrum: 

1 (PN dN (n-l)(n-2) 



2t: Pt dydpT dy 2TrnT(nT + ■mQ{n — 2)) 



m^-nia 



nT 



(14) 



The dashed-dotted curve in Fig. is the Levy function 
fit with ^ jndj — 0.90 to the K* spectrum in all the 
measured pt range (j)t < 4 GeV/c). 
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FIG. 11: (Color online) The invariant yields for both {K*° + 
'K^)/2 and {K*+ + K*-)/2 as a function of pr for \y\ < 
0.5 in minimum bias p + p interactions. The dotted curve is 
the fit to the power-law function from Equation 1131 for pr > 
0.5 GeV/c and extended to lower values of pr- The dashed 
curve is the K*^ spectrum fit to the exponential function from 
Eauation ll2l and extended to higher values of pr- The dashed- 
dotted curve is the fit to the Levy function from Eauation ll4l 
for Pt < 4 GeV/c. Errors are statistical only. 



C. Average Transverse Momentum (pt) 

In Au-|-Au collisions, the pt range of the exponential fit 
covers >85% of all the K* yield so that the K* average 
transverse momentum {{pt)) can be reasonably calcu- 
lated by using the inverse slope parameter (T) extracted 
from the exponential fit function and assuming the ex- 
ponential behavior over all the pt range: 



j^ple-^V^^o-"^-y^dpT 



In p + p collisions, the neutral and charged K* spec- 
trum shown in Fig.llllcovers >98% of all the K* yield so 
that the (jjt) is directly calculated from the data points 
in the spectrum. The systematic uncertainty in p + p 



includes the differences between this calculation and the 
exponential fit to the K*'^ only at pt < 1.6 GeV/c, the 
power-law fit to both neutral and charged K* at pT > 1-5 
GeV/c, the Levy function fit at pt < 4 GeV/c. The sys- 
tematic uncertainties for all the {pt) values include the 
effects discussed in the previous section and the differ- 
ences caused by different fit functions to the invariant 
yield, such as the Boltzmann fit (mTe"*^™^"™"^/^) and 
the blast wave model fit 0. The calculated K* {pt) for 
different centralities in Au-|-Au and minimum bias p + p 
collisions are listed in Tabic ITTTl 

TABLE III: The K* {pt) for different centralities in Au-|-Au 
and minimum bias p+p collisions. The first error is statistical, 
the second is systematic. 

(pt) (GcV/c) 



top 10% central 
0-10% 
10-30% 
30-50% 
50-80% 
p + p 



1.08±0.03±0.12 
1.08±0.08±0.12 
1.12±0.06±0.13 
1.08±0.05±0.12 
1.03±0.04±0.12 
0.81±0.02±0.14 
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FIG. 12: The K* (pt) as a function of dNch/dr] compared to 
that of 7r~ , , and p for minimum bias p-\-p (solid symbols) 
and Au-|-Au (open symbols) collisions. The errors shown are 
the quadratic sum of the statistical and systematic uncertain- 
ties. 

The K* {pt) as a function of the charged particle mul- 
tiplicity {dNch/dr]) is shown in Fig. 112! and compared to 
that of TT~ , K~, and p for different centralities in 
Au+Au and minimum bias p + p collisions. The K*'^ 
{pt) in Au-|-Au collisions is significantly larger than in 
minimum bias p + p collisions. No significant centrality 
dependence of {pt) is observed for K* in Au+Au colli- 
sions. 
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D. Particle Ratios 

The K* vector meson and its corresponding ground 
state, the K , have identical quark content in the context 
of the standard model of particles. They differ only in 
their masses and the relative orientation of their quark 
spins. Thus, the K* / K yield ratio may be the most inter- 
esting and the least model dependent ratio for studying 
the K* production properties and the freeze-out condi- 
tions in relativistic heavy- ion collisions. The K* and (j) 
mesons have a very small mass difference, their total spin 
difference is AS* — 0, and both are vector mesons. One 
significant difference between the K* and is their life- 
times, with the 4> meson lifetime being a factor of 10 
longer than that of the K* . Therefore, it is important to 
measure the (j)/ K* yield ratio and compare the potential 
differences in K* /K and (p/K yield ratios in relativistic 
heavy-ion collisions to study different hadronic interac- 
tion effects on different resonances. 



interactions at y/s^ = 200 GeV are compared to mea- 
surements in e+e" [48Ll49l [50 l l5l||.p-l-a [H^. and p -I- p 
[Hilil. The K*/K- yield ratios depicted in Fig. US) 
do not show a strong dependence on the colliding sys- 
tem or the cm. system energy, with the exception of the 
K* /K~ yield ratio at ^/s^ = 200 GeV. In this case, the 
K* /K" yield ratio for central Au+Au collisions is signif- 
icantly lower than the minimum bias p + p measurement 
at the same cm. system energy. The 4)/K* yield ratios 
as a function of the cm. system energies are depicted in 
the lower panel of Fig.^| The (p/K* yield ratios for cen- 
tral Au-f Au collisions at ^/s^ = 130 |3i| and 200 GeV 
and minimum bias p-\-p interactions at y's„„ = 200 GeV 
are compared to measurements in e+e" [s^, |^ 

and p + p m. ■ Figure El shows an increase of 
the yield ratio (t)/K* measured in Au-|-Au collisions com- 
pared to the measurements m. p + p and e"'"e~ at lower 
energies. 
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FIG. 13: The K*/K (upper panel) and 4>/K* (lower panel) 
yield ratios as a function of the cm. system energies. The 
yield ratios for central Au-(-Au collisions at -^/sTtn" ~ 130 |3fil | 
and 200 GeV and minimum bias p+p interactions at y^s^^iv ~ 
200 GeV are compared to measurements from e+e~ at -^i of 
10.45 GeV EIK 29 GeV and 91 GeV Ullsi, pp at -/s 
of 5.6 GeV lla and pp at y/s of 27.5 GeV 52.5 GeV 
m and 63 GeV The errors at y/sWN= 130 and 200 

GeV correspond to the quadratic sum of the statistical and 
systematic errors. 



The K* / K yield ratios as a function of the cm. sys- 
tem energies are shown in the upper panel of Fig. E| 
The K* /K~ yield ratios for central Au-|-Au collisions at 
yfs^ = 130 ,35.] and 200 GeV and minimum bias p + p 



TABLE IV: The K'/R-, 4>IK*, and (ji/R- yield ratios for 
different centralities in Au-|-Au and for minimum bias p + p 
interactions. The first error is statistical, the second is sys- 
tematic. 



K* IK 



0-5% 
0-10% 
10-30% 
30-50% 
50-80% 

P + P 



0.23±0.01±0.05 
0.24±0.02±0.05 
0.26±0.02±0.06 
0.26±0.02±0.05 
0.35±0.01±0.05 



0.60±0.06±0.12 
0.63±0.07±0.14 
0.58±0.06±0.13 
0.53±0.05±0.11 
0.53±0.03±0.09 



0.16±0.01±0.02 
0.15±0.01±0.02 
0.16±0.01±0.02 
0.16±0.01±0.02 
0.15±0.01±0.02 
0.14±0.01±0.02 



Table Irvl lists the K*/K-', (ji/K*, and (t^/R- yield 
ratios for different centralities in Au-t-Au and minimum 
bias p + p interactions. Figure H"!] depicts the K* /K~, 
't>/K~ ^^^d p^/tt^ 12] yield ratios as a function of 
dNch/drj at ^s„„ = 200 GeV. All yield ratios have been 
normalized to the corresponding yield ratio measured in 
minimum bias p + p collisions at the same ^s„„ and 
indicated by the solid line in Fig. 1141 As mentioned pre- 
viously and shown in Fig.^l the K*^ / K~ yield ratio for 
central Au-|-Au collisions is significantly lower than the 
minimum bias p + p measurement at the same cm. sys- 
tem energy. In addition, statistical model prediction of 
K*/K of 0.33±0.01 dlilil is considerably larger (in 
a 2cr effect) than than our measurement of 0.23±0.05 in 
0-10% Au-|-Au. The K*^ regeneration depends on ax-n 
while the rescattering of the daughter particles depends 
on (TTrTT and cr^jj, which are considerably larger (factor 
-5) than ax-n i22ti22,]. The lower K*°/K- yield ratio 
measured may be due to the rescattering of the K*^ de- 
cay products. The p^/tt" yield ratio from minimum bias 
p + p and peripheral Au-|-Au interactions at the same 
cm. system energy are comparable. Due to the rel- 
atively long lifetime of the (f) meson and the negligible 
axKi the rescattering of the decay products and the (p 
regeneration should be negligible. The statistical model 
calculations 0, predict the (/i/Tr" yield ratio to be 
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0.025±0.001 while STAR measured the K' /tt- yield ra- 
tio to be 0.15±0.02 [H. Thus the (ji/K' yield ratio 
combining the model prediction and experimental mea- 
surements is 0.17±0.02 which successfully reproduces the 
4>/K~ yield ratio measurement depicted in Table HVl and 

Fig. in 



E. Elliptic Anisotropy V2 

In non-central Au+Au collisions, the elliptic flow {112) 
is defined as the second harmonic coefficient of the 
Fourier expansion of the azimuthal particle distributions 
in momentum space ,531 • The K*^ V2 can be calculated 
as: 
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FIG. 14: The K* /K' , (j)/K' , and p° /tt' yield ratios as a 
function of dNch/drj for Au-|-Au collisions at y^s^Tw ~ 200 
GeV. All yield ratios have been normalized to the correspond- 
ing yield ratio measured in minimum bias p + p collisions at 
the same cm. system energy and indicated by the solid line. 
Both statistical and systematic uncertainties are shown. 



where (j) is the K*'^ azimuthal angle in the momentum 
space, ^'r- denotes the actual reaction plane angle and () 
indicates the average over all in all events. 

For each Ktt pair, the reaction plane angle was esti- 
mated by the event plane (^'2) which in turn was deter- 
mined by using all the primary tracks except the kaon 
and pion tracks in the pair: 

*2 = ^ X tan"^ 

J2i ^» sin(2(/>t) ~ UJK sm{2(f)K) - uJtv sin(20^) \ 
J2i cos(20i) - UJK cos(2(/)_ft:) - cos(20^) / ' 

where is the weight for each track used to optimize 
the event plane resolution, the subscripts K and tt stand 
for the kaon and pion candidate track, respectively. This 
prevents the auto correlation between the Kir azimuthal 
angle 4)k-k and the event plane angle ^2 113 • 



The centrality dependence of the resonance yield ra- 
tios depicted in Fig. 1141 suggests that the regeneration 
and the rescattering of the decay products are negli- 
gible, and the rescattering of the K*^ decay products is 
dominant over the K*'^ regeneration and therefore the re- 
action channel K* ^ Ktt is not in balance. As a result, 
the K*^/K~ yield ratio can be used to estimate the time 
between chemical and kinetic freeze-outs: 

I kinetic — j-^ | chemical X 6 , (f^) 

K K 

where r is the K* lifetime of 4 fm/c and At is the time be- 
tween chemical and kinetic freeze-outs. If we use the min- 
imum bias p + p measurement of the K*^ / K~ yield ratio 
as the one at chemical freeze-out and use the most cen- 
tral measurement of the K*^ /K~ yield ratio in Au-|-Au 
collisions for the production at kinetic freeze-out, then 
under the assumptions that i) all the K*s which decay 
before kinetic freeze-out are lost due to the rescattering 
effect and ii) there's no regeneration effect, the time be- 
tween chemical and kinetic freeze-outs is short and Ai — 
2 ± 1 fm/c. All the above assumptions reduce the esti- 
mated At. Thus the previous value is a lower limit of At 
and it is not in conflict with the estimations (>6 fm/c) 
in |42l |. These two measurements together indicate that 
considerable resonance regeneration effect may happen 
even (about 4 fm/c) after chemical freeze-out. 




Pt (GeV/c) 

FIG. 15: The K*° V2 (filled stars) as a function of pt for 
minimum bias Au-|-Au collisions compared to the (open 
triangles), A (open circles), and charged hadron (open dia- 
monds) V2- The errors shown are statistical only. 

In minimum bias Au-|-Au collisions, the unlike-sign 
and mixed-event Kit pair invariant mass distributions 
are reconstructed in cos[2((/) — ^'2)] bins and in pT bins. 
After the mixed-event background subtraction for each 
cos[2(0 — ^2)] bin and pt bin, the K*^ yields are then 
obtained as a function of cos [2(0 — \E'2)] for given px bin. 
The average (cos[2(0 — '^2)]) is then calculated for each 
Pt bin. The finite resolution of the event plane angle, 
which is due to the limited number of tracks in the event 
plane calculation, reduces the measured X*'^ V2. Thus 
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the above obtained {cos[2{(j) — ^2)]) values are further 
corrected for an event plane resolution factor (< 1) using 
the method presented in |^. Fig-urelT^shows the K*^ V2 
as a function of px compared to the Kg, A, and charged 
hadron V2 for minimum bias Au+Au collisions .27j ■ A 
significant non-zero K*'^ V2 is observed. Nevertheless, 
due to the large uncertainties on the K*^ V2 measure- 
ment, no significant difference is observed between the 
K*'^ V2 and the K'^, A, and charged hadron V2- 

In order to calculate the contributions to the K* pro- 
duction from either direct quark or hadron combinations, 
the following function SSj was used to fit the K*''' V2- 



V2{pT,ri) = 



1 + cxp[— — b)/c] 



dn. 



(18) 



where a, b, c and d are constants extracted by fitting to 
the Kg and A V2 data points in [s^. and n is the open 
parameter standing for the number of constituent quarks. 
From the fit to the K*'^ V2, n = 3 ± 2 was obtained. 
Due to the large statistical uncertainties, it is difficult to 
identify the K* production fractions from direct quark 
combinations {n = 2) or hadron combinations (n ~ 4). 
About 15-20 times more Au+Au collision events were 
taken by the STAR experiment in the fourth RHIC run 
in 2004 which is expected to provide enough sensitivity 
for more precise calculations of n to identify the K* from 
different production mechanism. 



F. Nuclear Modification Factor 

The number of binary collisions (Nbin) scaled central- 
ity ratio (Rcp) is a measure of the particle production 
dependence on the size and density of the collision sys- 
tem and is closely related to the nuclear modification 
factor {Raa)- Recent measurements of the A and Kg 
Rcp at RHIC have shown that in the intermedi- 



ate pt region (2 < pt < 4 GeV/c), the A and K^ Rcp 
are significantly smaller than unity. These measurements 
suggest that high pT jets lose energy through gluon radi- 
ation while traversing through dense matter. It has also 
been observed that the Rcp is significantly different for 
A and Kg with pr > 2 GeV/c. It is not clear whether 
this Rcp difference is due to a mass or a particle species 
effect. The K* is a meson but has a mass that is close 
to the A baryon mass. Thus, the measurement of the K* 
Rcp may help in discriminating between mass or particle 
species effect at the intermediate px region. 

The K* Rcp was obtained from the pr spectra of the 
top 10% and the 50-80% most peripheral Au-|-Au colli- 
sions. The K* Raa was calculated from the pr spectrum 
of the 10% most central Au-|-Au collisions and the pr 
spectrum of the minimum bias p + p collisions. 

The K* Raa and Rcp as a function of px compared 
to the A and if^ Rcp are shown in Fig. El The K* 
Raa and Rcp for pr < 1.6 GeV/c are smaller than the 
A and Kg Rcp indicating the strong rescattering of the 



K* daughters at low pT- The rescattering of the K* de- 
cay products is weaker for px > 1-6 GeV/c since K* with 
larger pT are more likely to decay outside the fireball [U . 
Therefore, larger pT K* have a larger probability to be 
measured compared to low pT K* . The K* Raa and 
Rcp are closer to the K^ Rcp and different from the A 
Rcp for pt > 1.6 GeV/c. Thus, a strong mass depen- 
dence of the nuclear modification factor is not supported 
and a baryon-meson effect is favored in the particle pro- 
duction in the intermediate pT region. 
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FIG. 16: The K* Raa (filled triangles) and Rcp (filled cir- 
cles) as a function of pr compared to the Kg (open circles) 
and A (open triangles) Rcp- The errors shown are statistical 
only. The dashed line represents the number of binary col- 
lisions scaling. The widths of the grey bands represent the 
systematic uncertainties of Raa (left) and Rcp (right) due 
to the model calculations of Nm,!- 



VI. CONCLUSION 

Results on the K*^ and K*^ resonance production in 
Au-|-Au and p + p collisions measured with the STAR ex- 
= 200 GeV were presented. The K*° 



KgTT^ at midra- 



periment at y/s^ 

and K*^ signals were reconstructed via their hadronic 
decay channels K*'^ Kn and K*^ 
pidity. 

The K*'^ /K yield ratios in Au+Au collisions were ob- 
served to be smaller than the ratio m p + p interactions 
which may be interpreted in the context of finite cross 
sections in a late hadronic phase. The result suggests 
that the rescattering of the K*''^ decay products is dom- 
inant over the K*'^ regeneration and therefore the re- 
action channel K* Kn is not in balance. As a re- 
sult, the K*'^ /K" yield ratio can be used to estimate 
the time between chemical and kinetic freeze-outs. Us- 
ing the K*^ /K~ yield ratio, the lower limit of the time 
between chemical and kinetic freeze-outs is estimated to 
be at least 2 ± 1 fm/c. 
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A significant non-zero K*^ elliptic flow V2 was mea- 
sured as a function of pt in minimum bias Au-|- Au colli- 
sions. Due to limited statistics, no conclusive statement 
can be made about the difference between the K*'^ V2 and 
the Kg, A, and charged hadron V2. The estimated num- 
ber of constituent quarks for the K*'^ from the V2 scaling 
according to Eauationll8lis 3 ± 2. Thus, larger statistics 
for Au-t-Au collision data are needed to identify the K* 
production fractions from direct quark combinations or 
hadron combinations. 

The nuclear modification factors Raa and Rcp 
were measured as a function of pr- Both the K*^ Raa 
and Rcp are found to be closer to the Kg Rcp and 
different from the A Rcp for pr > 2 GeV/c. A strong 
mass dependence of the nuclear modification factor is not 
observed. This establishes a baryon-meson effect over a 
mass effect in the particle production at the intermediate 
Pt region. 
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